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NATIOHAL AIJVTSCKy COMMIOTEE FOE AERONAUTICS 


TECHNICAL NOTE NO- 1381(- 


A REVIEW OE BOUNBARY-LAYER LITERATURE 
By Neal Tetervin 


SU14IAEY 


A concise noranatheaiatlcal review of the subject of boundar,y 
layers is presented. The contents, althou^ insufficient foz’ the 
solution of specific problems, are sufficient fcr an introduction 
to the subject. A. list of reference papers is given frcan. -vdiich the 
detailed Isnowledge necessary for the solution of specific problems 
can be obtained. 


INTEOLTJCTION 

The literature on the subject of bouru?J3ry layers contains so 
many papers of varying quality that it is difficult for a- newcomer 
to the subject to choose the papers that provide the maximum gain in 
knowledge for the effort expended. Althou^ references 1 and 2 pro- 
vide detailed reviews of bomdary-layer theory, no short noranathe- 
matical summary is readily available. 

The purpose of the present paper is to provide a short sunmiary 
that contains exact or approximate information that is believed to 
be i^seful. The summary is confined to cases for -v&ilob. the physical 
properties of the fluid are constant, that is, to incQmpres'Bible 
flow with no temperatiJre effects • An introduction to boundtoy-layer 
literature is provided, and reference papers are listed from, which 
information on subjects of special interest may be obtained. 

The material presented herein was originally presented os a 
talk at Wrlfjht Field, Dayton, Ohio, on September 26, 19^5 . 

BOUNDARY-LAYER TEEMINOLO&Y 
Symbols 


p density 

u velocity couponent parallel to surface 

V velocity ccmponent porpondicular to surface 
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S 

X 

y 

p 

1^ 

TJ 

S 


8* 

8 

H 

at 

% 

{52 

{ dx 


distance along sitrface from leading edge 
distance normal to surface 
statlo lOTessure 
coefficient of viscosity 

velocity can3)onent parallel to surface at outer edge of 
"boundary layer 


momeritum thicloaess 



displacement thiclsness 



nominal thiclcnese of "boundary layer 
surface shearing stress 


ratio of displacement thlcloiess to momentum thickness 



dynamic pressiire 



fe)(?r 


kinematic viscosity 
free'-stream velocity 

absolute magiltude of rate of change of U irtth x 
Reynolds number 
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T 

I 


chord 

Reynolds nimifaer 


(¥) 


local, shearing stress 
mlzing length 




m 

esponent 


value of 

T!i= 1 

-iT) 


Reynolds 

n 

eaponent 


/USX 


eajjonent in formula for variation of velocity along surface 
dx 

V _ X dU 

■ ■ 

k constant in formula for variation of velocity along surface 

nt 





u'v' dt 


VO 


u' x-component of fluctuation velocity 
v' y-component of fluctuation velocity 
t time 

K von Ka3r233an's universal constant 
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D 

6 

X 

c 


Eeynolds number 

height of r outness per tide 
angular velocity 
radius of disk 

thicloiess of laminar sublayer 


constant in surface-fi-iction formula 
constant in surface -friction form'da 


Pohlhausen shape parameter 



inflow velocity 


Terminology 

A boundary layer may be defined as a region in the flow field 
in Tdiich the viscous forces in 'bhe equation of moticai are not all 
negligible. (See appendix.) For flow over bodies at the Reynolds 
numbers encountered in applied aerodynamics^ the region in which the 
viscous forces are not ne^igible is confined to a thin layer of 
fluid next to the surface, the Prandtl boundary layer (reference 1) . 
The viscous forces are confined to this thin boundary layer because 
the space rates of change of shearing stress are largo enough to 
produce other than negligible viscous forces only^ in the thin surface 
layer of fluid. In the absence of solid boundaries, boundary layers 
occur where streams of fluid that move with different velocities are 
in contact^ for 0 xan 5 )le, Jets and wakes. Boundary layers may be 
divided into two classes, laird, nar and turbtilent. • 

A laminar boundary layer is one in which the paths of the 
particles of fluid never cross one anotherj the nei^boring layers 
of fluid glide over one another as if they were solid sheets and all 
interchange of momentum between adjacent layers takes place only by 
molecular motions. A turbulent boundary layer, on the othor hand, 
is on© in which the paths of the particles of fluid- cross one another 
and in which almost all of the interchange of momentum between 
adjacent layers is caused by the irregular motion of small fluid 
masses. 
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Tn order lio dlocuss the boundaa:^ layer It Is first necessary to 
define the terras velocity profile^ surface fh’lction./ boun'iary -layer 
thiclaiess, and separation point. A velocity profile is the curve 
that gives the distribution of the component of velocity parallel to 
the surface "jd-th distance natrmal to the surface (fi£« 1) . In bovndarj’’ 
layer theory th© velocity ccaaponent parallel to the surface is ©q.ual 
to the ma^itude of the total velocity because the velocity component 
norroal to the surface is negligible. The surfs c© friction is the 
shearing stress between the fluid and the solid body. 

The thickness of the boundary layer may be defined as the dis- 
tance^ in direction normal to, the eurfacey at which the total pres- 
sure (■•f) differs by an arbitrary small amount from the total 

pressure of the undisturbed flow# At distances from the body greater 
than S, the flow Is assimied to be inviscid (fig. 2) . 

The separation point is the point on the siTfaco of the body at 
tdiich the surface friction is zei'o. ppstream of the point the direc- 
tion of flow In the boimdary layer next to the s'.xface is dotejatream, 
and downstream of the point the direction of flaw ?n the boundary 
layer next to the siarface is upstream (fig. 3)< 


LAMIHAR BOUiroAEY TA'dSP. 
Methods of Calculation 


By the use of the boundary -layer og,uations of motion together 
vrLth the conditions that the solutions of the eq.uationo mvist satisfy 
at the innex’ and outer edges of the boundary layer tlie character- 
istics of the laminar boundary layer over a bodj'- may be determined 
completely when the velocity distribution over the body outside the 
boundary layer is known and the flow is such that ’ulie boundary-3nyer 
approximations are applicable. In order to avoid the purely mathe- 
matical difficulties associated nrilth the exact method of solution, 
various approximate methods for the ocaneputation of the velocity 
profile, surface friction, boundary-layer thickness, and separation 
point have been developed. 

^ f ft 

Yon Karman mcmentum eq.uation «- The von. Earman moment vm eq.ua- 
tion (reference 3) results ftrom the application of. the momentimi 
theorem to boundary -layer flo\r and mahes it possible to compute the 
boundary-layer thickness over a body ■jrf.th an arbitrary pressure dis- 
tribution -viiethor the flow in the boundary layer is laminar or 
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turbuleht if the(Taridtlon of the ■boimdai:‘y_‘^^yer velocity profile 
and surface friction over the "body is knovttxi 

the eq.uation is obtainable by tw inethods; by applying the 
momentum principle to a box^ the lo'tver side of Tillich extends for an 
infinitesimal dis-bance alcaig -the solid surface, the upper side. of 
Tdiich is the nominal thickness of the boundary layer, end the sides of 
Tdiioh are planes perpendicular to the body surface (appendix), or 
by integrating the boundary^layer equation of motion ^-with respect 
to -the distance normal to the surface. The von Karmen momentum 
equation contains the same assumptions as the Prandtl boundary-layer 
equations (appendix) and may be liritten as (appendix) 


p ^ = s* ^ + T 

dx dx ° 


( 1 ) 


The quantity 0, Tdiere 


pt ■ pd 

pU ^0 s pU I u dy - p 

Jo Jo 


dy 


is called the momentum -thickness and is a length of such magaitude 

that represents the difference betTreen the rate of momentum 

flow that would exist if the mass floT-dng -throu^i a boundary-layer 
cross section liere moving with the velocity at the boundary-layer 
edge and the actual rate of momentum flow tlirou^ -the bouiidary-layer 
cross section. 

The quantity 6*, uhero 


P5 


pus* = pus - p 


uo 


u dy 


is called the displacement thickness and is a length of such magni- 
tude that fluid floTdng tiirou^ it Td.th the boundary-layer edge 
velocity U produces a rate of mass flow equal to -the difference 
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between the rate of mass flow which would exist if the fltiid through 
a cross section of the boundary layer were flowing with the houndary- 
layer edge velocity and the actual rate of mass flow through the 
cross section of the boundary layer. 

Eq;uatlon (1) indicates that the effect of pressure gradient on 
the momentum, defect of the fluid flotfing through the boundary layer 
is directly proportional to the displacement thickness of the boundary 
layer . 

For use in the Pohlhausen method and for the coaiputation of 
boundary-layer thicktiesses; the von Karman equation is written as 


dx 



U dx 


( 2 ) 


Pohlhausen *s method .- The Pohlhatisen method,, an approximate 
method baaed on the von Karman momentum equation has been vriLdely 
used and is useful for obtaining qualitative information. The 
purpose of the Pohlausen method (reference 1, pp. 103-ll g ) is to 
ccaaput© the characteristics of the laminar boundary layer in two- 
dimensional flow when the pressure distribution outside the boundary 
layer is a known function of the distance along the surface. The 
method is based on the assumption that all laminar boundary-layer 
velocity profiles are given by a fourth -degree polynomial. Pohlhausen 
chose a fourth-degree polynomla3. after trying first-, second-, and 
third-degree polynomials, because the fourth-degree polynomial gave 
better agreement between his method and the exact Blasius solution 
for the velocity profile and skin friction on a flat plate than 
polynomials of lower degree. The coefficients of the fourth-degree 
polynomial are chosen to make the equation for the velocity profile 
satisfy the boundary-layer equation *f moti*«i at the inner and outer 
edges of the boundary layer. The result. is an equation for the 


velocity profile ^ = f 


•mere X = 


(See appendix.) 


u dx' 

The velocity profiles are thus a single -parameter famJ.ly of cvrvea 
in which the parameter X depends on the previous history of the 
botindary layer only through the thickness 6. The paraiaeter X is 

dx 


proportional to 


(appendix), the ratio of the unbalanced 
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horizontal presaxnre force acting on the small hox nsed in deriving 
the von Kannan nouentutix equation to the shearing stress acting on 
the em’face side of the box. The parameter X can be shown to be 
independent of (appendix) . For the interior of the boundary 

layer the equation of motion is ignored, but the von Kamaii momentum 
equation is satisfied, By satisfying the von Karman momentum equa- 
tion, the distribution over the body surface of tlie velocity profile, 
boundary- layer,, thickhes 8, . surface friction, and surface -pressure 
distribution is made consistent with the momentum theorem. 

To obtain the boundary -layer characteristics over the surface, 
the differential equation resulting from the substitution of t!ie 
equation for the velocitj’" profile into the von Karman aquation is 
solved (reference 1, pp. IO8-II2) . The method does not suffer 
from serious InaccuracioB for 0 < X < 12. For flow over a flat 
plate, X = 0, the skin friction differs from the exact value by 

onlj'- 3~ percent. Vhen, however, the pressure outside the boundary 

layer rises in the direction of floif, X < 0, the method becemes 
inaccurate. An investigation of -tii© reason for the inaccviracy of 
the Pohlhausen method (reference 4) leads to the conclusion that an 
inherent characteristic of the method is the late prediction of the 
separation point because the fourth-degree poiynamial_for the velocity 
distribution is not a good meth^Batical substitute for- the actual . 
velocity profiles obtained in the exact solutions. Althcuj^ the 
method is inaccurate in an adverse pressure gradient, it may often 
be used to obtain rapidly q-ualitative results concerning tiie effects 
on the velocity profile and surface friction of. changes in the pres- 
sure dlsti'ibution or boimdary-layer thickness (appendix). The method 
is an example of the fact that a theory wtiich ignores the. equation of 
motion in the interior of the boundary layer and therefore neglects 
the acceleration terms in the equation of motion will make the 
boundary -layer profile dependent on only the local conditions . A 
coaqjlete theory would make the space rate of change of boundary - 
layer profile, rather than the profile itself, depend on the local 
conditions . 

Falkner 's method. . - The purpose of Falkner 's method (references k 
and 5) is to provide for the practical coinputation of the surface 
friction, momentum thlcknese, and displacement thickness of .any two- 
dimensional laminar boundary layer. The method is based on existing 
tables of- solutions of the boundary -layer equation for special types 
of pressure distribution (reference 4) . Xn order to obtain the 

tables, Falkner expanded ttie quantity 
series from the stagnation point u’itla 


■( 


-Vaf* 


' o 

pjj2j \-o J 


in a Taylor 
X as the Independant variable . 
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The coefficients in the Taylor series are the deriyativee of ^ 
i/ith respect to x and were obtained from the derivatives of the 
complete houndary-layer eq.uatioh at the stagnation point, with the 
help of Isnowi solutions of the special eq.uation (reference 4) to 
which the ccmplete eq.uation reduces -in the vicinity of the stagna- 
tion point. The series espansion for ^ tras then idicited to a 

'*1 ^ 

special type of pressure distribution given by TJ = Hod o , 
where X’ is a constanb. The separation point = 0) for the 
special type of pressure distribution was determl.ned by using the 
Taylor series expansions for After the separation points were 

computed, the solutions for different pressure distributions of the 
same family were tabulated and completed^ these tables form 
basis of the simplified laethod of calculation (reference 5 ) • 

The assumptions in Fa Umar's simplified meiiiod are: (l) The 

special forjta of the boundary -layer equation of motion represents 
the conditions of flow near the stagnation point with good accuracy. 
(2) The surface friction at a point is given accurately by replacing 
the actual pressure distribution by one of tiie particular famj.ly of 
pressure distributions; a new px’essure distribution is chosen for 

each point. ( 3 ) The relation H = — is a function only of 

Q 

and this function can bo determined from solutions of the boundary- 
layer equation of motion near the stagnation point. 

The computation that must be made to solve a problem is simple 
and rapid. It consists of evaluating a simple integral (reference 5) 
from the given data and using the results of the integrations with 
the standard tables to obtain all the q.uantltiea that are of interest 

The method is shown to be sviitable for the computation of the 
separation point by the good agreement between the ccmqputed separa- 
tion point and the known separation point for two cases: One, an 

exact solution of the bc\jndary-layer eq.uations and the other, an 
experimental determination of the separation point (reference 5 ) • 

Hartree ' s method . - The purpose of Hartree's method is to pro- 
vide an accurate method fear the computation of all the di^acter- 
istics of the laminar boundary layer. The basis of the .met-hod is 
the integration of the boundary -layer, equation of motion by the 
replacement of tlie equation of motion, a partial differential equa- 
tion, by an equation involving finite differences and ordinary 
derivatives. The only approximation used, other than those con- 
tained in the boundary-layer equation, is tloat a derivative may be 
replaced by a finite difference. 
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To order to compute the 'boundary-layer characteristics from the 
given pressure distrihutlon, the partial differential e(j.uation of 
motion is replaced by an approximately equivalent ordln^y differ- 
ential equation by replacing the derivatives with respect to one of 
the variables by corresponding finite-difference, ratios. The derlva 
tives with respect to the other variable are integrated either 
mechanically or by some standard process for the numerical inte- 
gration of ordinary differential oqijations. Derivatives paraJJLel 
to the boundary are replaced by finite differences, and integratlcai, 
is carried out along successive normals to the boundary at finj.te 
Intervals so that from the distribution of velocity across one 
section of the boundary layer the distribution of velocity across 
another section at an interval downstream is calculated. The limit 
to the accuracy obtainable with the method is the amount of work, 
which increases with the accuracy desired. 


Critical Eemarlcs 

The Pohlhausen method is useful f or obtaining goneral_gTjall- 
tative infbi^^lm for~*e'lth er "ra?orable~ or adve rse p ressure gradients. 
and'fbr ^obtaining quantitativ e informat ion ' for ifayoreble pressure . 
g'a dlents ^ falkner 's method seems to provide sufficient accuracy 
for the solution of problems and to be rapid. Eartree 's met^iod is 
potentially more Impoirbant than any of the others because if the 
ccmrputational work can be decreased the method can provide an 
accurate solution of the bovindary -layer equations. At present .the 
method is useful for providing solutions for testing apisroxlmate 
methods. The methods of references 6 to 12 seem to be inferior to 
Falkner 's rapid method for general vise because they either take 
much longer or are not -so accurate. 

■When rapidity of computation is of prime Irqportance, tlie posi- 
tion of the separation point of the laminar, boundary layer may be 
estimated by replacing the velocity distribution over the body In 
the region of adverse gradient by a velocity distribution with a 
constant gradient. Approximate methods besed on the use of a 
constant gradient may be obtained from the Polxlhausen method (refer- 
ence 1, pp. 108-112), frcm the von Karman-Mi.llikan method (refer- 
ences 7 and 8) , and from tha work of Howarth (reference 6) . A 
method based on Pohlhausen 's work will usually predict separation 
too far downstream, and the methods based on the von Kaiman-Mlllikan 
method will usually predict separation too far upstream. The 
accuracy to be expected frcan a method based on Hovrarth's exact 
solution is unknown but should be good when the actual gradient 
is close to constant. 
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Exact Solutions of Boundary-Iiavei* Eo nations 

Biasing solution of 'bound.aiy -layer equation for flow over a 
flat plate -with zero pressure gradient .- The results obtained fi’om 
the Blasius solution (reference 1, pp. Sk^go) are 

( 1 ) Th e thiclmess of the boundary layer Is proportional 
to (/xP/Uq . 

(2) The dreig coefficient based on the dra^ of one side of a 

pL 

flat plate and the projected area of the plate is Cjj = — ~r-» 

(3) -The velocity ratio u/TJq in the boiandary layer is a 

fimction of the single variable ti = . — , 

V^Wo 

Hovarth .- Howartli (reference o) solved the boundary- layer 
equations of motion for the case of flow ever a flat plate with the 
velocity outside the boundary layer decreasing linearly in the dlrec 
tion of flow and with zero thichness of the boundary layer at the 
plate leading edge (figf if-). The restilts obtained were 

( 1 ) The thickness of the boiindary layer depends only on 
and on j2|g-. 

( 2 ) The velocity ratio u/Ug in the boundary layer depends 

only on — 7— - and on {— j 

l^zu/Uo I<i3c 1 Uq 

(3) The local surface-fricticai coefficient To/Sq decreased 
from an extremely large value at the leading edge to zero at the 
separation point. 

(h) The amount of velocity recovery AU/u^j before separation 
(fig. 4-) is a constant and is independent of the Reynolds number end 
of the rapidity with which the velocity is recovered. The rapidity 
of velocity recovery does not appear because the initial boundary- 
layer thickness is zero. 
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Yon E^man .- Yon TSsrwan reduced the con^leto equations of 
motion for a rotating disk in laminar flow to a system of ordinary 
differential -equations and obtained approximate -solutions (refer- 
ence 3) • The results 'tfere 


(l) The boundary “layer thickness is constant over the rotating 


disk and is 



(2) The inflow velocity normal to and far from the disk 
is C = 0 .J0Q\j\xa. 

(3) The turning moment required to rotate the disk for both 

sides is M = a^p 

. coa^ ^ ^ 

V u 


Cochran (reference 2) solved the system of ordinary' differential 
equations exactly by a numerical process. As in any exact solution 
of the equations of motion or of the equations of the boundary layer, 
no. definite bomdary-layer thickness was obtained. The constaht in 
result (2) for C was found to be 0.866 instead of O.708 and the 
constant in result (3) for M was found to be 3.67 instead of 3.68. 


Conclusions from Frandtl boundary-layer equation .- Useful 
information can be obtained from the Prandtl boundary -layer equa- 
tion as given by Falkner in reference 4 without obtaining solutions 
of the equation. The conclusions are 

( 1 ) For a fixed velocity distribution along the body and- a 
fixed point on the- body, the nondlmensional thickness 6/c of the 

boundary layer is inversely proportional to 

( 2 ) For. a fixed velocity distribution along '.the body and a 

fixed point on the body, the local starface-frlction coefficient — 

ju~C r— ' 

is inversely proportional to u-^ = total drag coeffi- 

cient of the part of a body covered by a laminar boundary layer 
therefore varies as 1 /Jr^. 


TJ^ 

u 


(fig. 5) V 



NACA TH No. 138^1- 


13 


(3) J'or a given velocity distribution over tiie body, the 
separation point is independent of the Reynolds number Rq- This 

fact is valuable in esperliaental irork £ind in cases -where coiaputa- 
tions of boundary -layer cliaraoteristics are made for mr.-re than one 
Reynolds number . 

(if) For a fixed velocity dia-bribution along the body and a 
fixed point on the body, the curve of vl/U against y/S is 
invariable and is independent of -the Reynolds number. The curve 

of u/tr against — \Jr^ is invariable ; this provides a good method 

for testing vhether a velocity distribution is' laminar. 


TURBULEHT BOUITOARY LAYERS 


In contrast to the smooth flow associated -wi-fch laminar boundary 
layers, a mixing flow is generally associa-fced -^dLth turbulent boundary 
layers. In a turbulent boundary layer the momentxm interchange 
bet-ween adjacent fluid layers is caused mainl^; by the irregular 
motion of small fluid masses'. The information concerning turbulent 
flow is largely empirical; -vhereas the information fo^- laminar motion 
is obtained -wholly from tho equations of motion. 

Skin friction .- Empirical skin-friction form-ulas for flow over 
smooth flat plates -wl-bh zero pressure gradient. and for flow in pipes 
-where a small favorable pressure gradient exists are available 
(reference 1 , pp. 135 - 15 il, and references I3 to 15) • Empirical 
skin-friction formulas are also given by Goldstein in a British 
paper of limited distribution. These references s^iow that 

(1) For eq.vtal Reynolds numbers the turbilLent skin-friction 
coefficient is greater than the laminar skln-fr-’.ctlon coefficient. 

(2) The tvirbulent skin-friction coefficient decreases less 
rapidly than the laminar skin-friction coefficient as the Reynolds 
number increases . 


The tixrbulent skin-fric-bion coefficient can be calculated from 


the same relation as for laminar flow 



if the velocity 


profile is known inside the laminar sublayer, that is, the region at the 
™il Id. -which the velocity- fluctuations disappear and the flow is laminar 
(reference 18) . An estimate of the thickness of -the laminar sublayez’, 
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■based on data from flow tiircuf^ pipes, la t^iven bj the relation 

~ ■ (reference 15 ) « 

Shearing stress »■ No exact relation Is kno\jn between the deriva- 
tive of the average velocity at a point inside a turbul.ent boui^ry 
layer and the shearirig stress at the .point. The absence of this 
inf carnation forces recourse to experiment to obtain iaiformation con- 
cerning turbulent flow. If a relation betw-een the local shearing 
stress and the derivative of the average local velocity is assumed 
a method for computing the diaracteristics of the turbulent _bo’u.dary 
layer may be devised. The method, however, will depend on the 
assumption concerning the sheering stress. 

Prcm the fun damen tal relation for the shearing stress in turbu- 
lent flow T = -pu'v' (reference 1, pp. 119~134, and reference 15), 
Prandtl derived the mixing-length theory, an approximate theory that 
relates the local shearing stress to the local density, the local 
derivative of the average velocity, and tlio local value of a so-called 

mixing length I (reference I 5 ) . Th® equation is t = pZ j 

\Sy/ 

The mixing-length theory is baaed on the fact that if two adjacent 
layers of fliiid have different velocities, the interchange between 
the two layers of small masses of fluid that have the velocities of 
the layers frcaa which they come will tend to eq.uallze the velocities 
of the fluid layers. The tendency toward equalization of veloclti’es 
may be considered as being caused by an apparent shearing stress 
between the tw layers. Yon Karinan (reference 15) obtained an 
approximate expression for the mixing length as a function of the 
local velocity derivatives and a universal constant by assuming 
that the la^ocess of turbulence at ary two points Is similar and 
differs only in the length and- time scales. The relation can be 
used to compute a velocity distribution in pipes that agrees very 
well with the experimental velocity distribution. 

Boimdary -layer velocity profile '. - The velocity profile for the 
turbulent boundary layer differs markedly in appearance from the 
velocity pc*ofile of the laminar borindary laj'er (fig. 6).. At large 
boundary-layer Eeynolde numbers the turbulent velocity profile shows 
an extremely rapid rise- in velocity in a very short distance. This 
lar^ slope at the wall causes the turbulent skin-friction coeffi- 
cient to be higher than the laminar skin-friction coefficient. 

■When the pressure gradient along the surface i.s zero or very 
subtly favorable, such as on plates or in pipes, the velocity 
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profile can "be approximated bj a logarithmic curve. Tlie ©ouation 
is (reference 15 ) 


= A log ^ + B 
V . t) 


and follows from, the assuaptiona that the shearing stress is constant 
across the pipe or boundary layer and eq.vial to the wall shearing 
stress and that the mixing length is given by the pq.uation 

I = Ey- 

For a wider range of pressure gradients including those which are 
adverse (that is, those in which the static pressure rises in the 
direction of flow), the velocity profile can be approximated by a 
power cvtrve 


u _ /y 

u \ q ) 


Both of these eq.uations for the velocity distribution become Inaccu" 
rate at distances from the wall ccmparable to the thichness of- the 
ifl’minRT ’ sublayer. At the wall both eq.uatlons Incorrectly give 
infinite values of du/dy. 


If the flow continues in an adverse pressure gradient, the 
velocity profile undergoes a change from one having hi^ velocities 
near the surface to one having low velocities near the surface 
(fig. 7 and reference 17 ) • 


For some purposes it seems permissible to assume that' the 
turbulent boundary -layer profiles form a single -parameter family 

5* , 

of curves and to • us© the ratio H = — , or the velocity ratio u/U 

e 


at some fraction of the houndary-layer thickness from the surface, 
as the parameter. The turbulent -velocity profiles usually found in 
xindisturb^ flow are simple curves and may be specified by . any one 
of a number of sui-babl;v chosen parameters. The suggestion that 
turbulent boundary- layer profiles form a single -parameter family of 
curves first appeared in a paper by G-ruschwitz (reference 18) in 


■vdiich the factor 



was used as a parameter. 
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Separation ^olnt ♦ - In an adverse pressm’e ©radjont the turbu- 
lent boundary layer will eventually separate* Por similar condi- 
tions, however, the percentage of the initial dynamic pressm-e 
wiilch can be converted into static pressure in flow over ettrfacos 
covered by turbulent boiandary layers is greater than in flow over 
surfaces covered by laminar boundary layers. Because the mixing 
between the inner and outer layers of fluid is much greater in a 
turbulent boundarj'- .layer than in a laminar boundar;/- layer,- more 
downs-broam momentm is transferred .from the outer layers of fi.uid 
to the inner layers of fltiid. 5 ?he towns bream mcmentura of the 
fluid near the wall in turbulent flew is -therefore maintained in 
easea_ln.’[dilch the doimst-ream mnm^antum of laminar fl ow woiU.d be 
e:idiaueted by the surface shear and separation would occur. 

Only empirical methods are available for estimating the 
separation point of turbulent boundary i^'ers. Of all the methods 
(reference , 1, pp. 155 '1^2; references it to and Garner's meth<3d 
which is given in a British paper of limited distribution) the two 
which seem to be most useful for pi'edictlng the behavior of the 
turbulent boundary- layer are those of Garner and reference I 7 . T-he 
Gruschwltz method (reference 18 )' introduced Jhe idea of a sin^.e- 
parameter family of ciirves for the velocity profiles of the turbu- 
lent boxmdary layer and made the rate of change alca^g the surface 
of the shape parameter rather than the shape parameter itself depend 
on the local '■conditions but did not make possible the ccanpu-batlcn 
of the separation point with sufficient accuracy for engine ering 
use (reference 23) . An attempted improvement o-f the- Grixsch'^-dtz 
method by Eehl (reference 20) has not been tested for- its ability 
to predict the separation point. 

The method of reference I 7 seems to be the most reliable method 
available at present for -the estimation of -bhe separation point of 
the turbulent boundary layer. The method uses two eq.ua-blDnB to- 
deteiTBine the behavior of a turbtilent boundary layer. The first 
equation is the von Karman momentum equation (equation (2)) . The 
second equation is an empirical equation -that gives the rate of 
change of boundary -layer shape parameter along the surface as a 
fmetion of the local conditions. The equation for the rate. of 
change of boundary-layer shape parameter was dovelcpod in the 
following manner: 

It was first yei'ifiod -that fer tlie eaperlr'.e:/;;^! data avar'.lable, 
the velocity profiles of the tvsrbulent boundary layiK.' formed a. 
single -^parameter family of curves idth H “■f(qq) as the .parameter. 

The assumption ■was -(dien-made that -the rate of change of boundary- 
layer shape parameter is a functicsn of the ratio of -the local 
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pressure gradient — to.ijb© local s]sin:’-friction coefficient t J2q, 

a dz" £ M ■ *■ 

and also to the local value of H. The ratio ■2;-'^ has the same 


physical significance as the parameter X in the Poh^ausen a»thod*' 

The q.uantities ^ ^ . and H were determined front the available 
q. dx 

experimental dataj the term ’ calculated the.sfcir- 

friction formula of reference 24. • From analj'^sis of thh extJerimental 

<3^ Q do ^0 » - ' ■ 

data, the variation of 9 — with — ~ and E can be repre- 

02 q. dx Tq 

sented by the eg.uation . 


. g ^ ^4.68o(E'’2,975) 

dx 

Equation (3) is solved simultaneousiy with the von Karman mcmentuia 
equation by numerical methods. For a fixed pressure distribution 
and transition point the method indicates a sllf^t forward movement • 
of the separation point with increase in Reynolds number. 

The method of Gamer differs from that of reference 1? orf.y by 
the use of a different empirical skin-friction relation and by tlie 
use of different constants in equation (3)* T!he different constants 
were obtained by analyzing the experimental data in reference 17 
and adding a small amount of data from experlmonts by Burl and 
Nikuradse. 

None of the msthods give any informatian on the variation of 
surface friction with boundary -layer prof ile ’ shape j all use' empirical 
skin-frlctioTi formulas derived from experiments with flat plates. 

Momentum thicknes s.- An approximate method fcsr tlie computation 
of boundary -layer momentum thicknesses that is useful ffjr the esti- 
mation of full boundary-layer thldoiesses arid profile -drag coeffi- 
cients when flow separation is not involved is given in reference 25* 
The method is based on the fact that if H is fixed at an average 

■ T” Q ^ 

value and that if a skin-friction equation of the form j-^.= 

IS used, the von KarHsn momentimi equation can be integrated.. The 
resTolt is a ormula for the computation of 9 . in flews with pres - 
sure gradients. 


- 2.035 (H - 1.286) (3) 

4 dx To . _ 
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The Bame type of forniijila is applicable to laminar boundary 
layers ■when the average velocity gradient is small. The fonrjula 
for laminar boundary layers has also been given in reference 26. 


Roughness . - Empirical skin-friction formuLae for flov over 
roug^i flat' plates and for flow in rou^ pipes are available for 
certain types of roughness (reference 1, pp. l45“154 and refer- 
ences 2, 15, and 27) . A characteristic of the flow over rough 
surfaces -with zero or small pressure gradient is that beyond a 
certain Reynolds number, •which depend-S 'csi the roughness, the^ skin- 
friction coefficient becomes Independent of the Reynolds number 
(reference 28) . The skin-friction coefficient of rou^ rotating 
cylinders (reference 3^) 'wa® found to become cons'bant at suffi- 
ciently large Reynolds nuiabers for a saturation density of roughxiess 
particles. Eor other than saturation densities, -the drag coefficie:it 
was concluded to decrease '^■dth Reynolds number. 


The addition of roug^moss to a surface covered by a turbulent 
boundary layer increases •the drag coefficient ■tdion the roughness 
hei^t becomes ccaaparable with the height of the laminar sublayer. 
The addition of roughness to a smooth plate will have no effect on 


the surface friction if 



•vhere 


,/s» 


is the surface" 


friction coefficient for the smooth plate (reference 30) 


TRANSITION 


A body in a stream usually has a laminar bcalndary layer for 
some distance from the stagnation point and behind the laminar 
boundary layer, a turbulent-boundary layer that extends to the 
trailing edge. The process of change of the laminar boundary layer 
to the tijrbulent boundary layer is fcnotra as transition. The appear- 
ance of the turbulent type of flow can be detected by the appearance 
of random fluctuations in the velocity components, by the change in 
velocity profile from one having a gradual increase in velocity 
■with distance from the •walls to one having a. much more rapid rise 
in velocity (fig. 6), and by the increase in skin-i’rietion coeffi- 
cient} a grea-fcer increase occurs in skin-friction coefficient at 
large Reynolds numbers -than at small Reynolds numbers . 

The flo’w conditions -that are favorable for the delay of -bransl- 
tion are: a , small Reynolds number, freedcsm from disturbances, and 

static pressure decreasing In the direction of flow (references 31 
and 32) . An investigation of the effect of curvature (reference 33 ) 
results in the following conclusions: The mechanism of transition 

is different on concave and convex walls, ■the transition poiii-fe^s 
not affected by convex curvature, concave ourvatui’e has a s-brong 
dos-babildzing effect, ■the influence of pressure gradient on transi- 
tion on a concave ■wall is negligible, the effect of pressvire gradient 
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on transition on convez rails is strongest near zero, pressure 
gradient, and stream tarbulenco has atout the sane effect on ti'unsi- 
tion for flow over concave and convex rails. 

■JOien the static pressure increases in the direction of flow, 
the maxi 3 p.um poseihlo length of laminar flow is- thy di.atance "between 
tile minimum pressure point and the Itiminar separation point. "Whether 
transition or laminar separation occurs first clopende on the Eejoaolds 
numTcer, the disturhances to the flow, and the strength of the adverse 
pressure gradient. As the adverse pressure gradient becomes smaller, 
the livelihood of transition occurring before sopareticn becomes 
greater; the flat plate represents the extreme case in which separa- 
tion never occurs. 

The instability of fee laminar boundarj'" layer on a flat plate 
has been investigated theoretically (refci'ence 3^) and the essentials 
of the theory have been verified eaperlmentally (reference 35) • A 
process for the computation of the instability point of the laminar 
boundary layer in the presence of pressure gradients, with examples 
of the results for airfoils, is given in reference 36 . Because it 
takes some distance for the laminar flow to become turbulent after ■ 
passing the instability point, the trarisition point is doxvnstream 
of the instability point. A congiariBon between the expex’imental 
transition points and the theoretical instability points for an 
WACA airfoil is given in references 36 and 37 • The ccmiparlson shows 
that, as is esepeoted from the theory, a deci’easing static, pressure 
in the direction of flow causes the distance bet^vroen the instability 
and transition points to increase and also, causes the instability 
point t(> move farther downstream at a fixed Eeynolds number. 


EEG-IONS OF SSPARATED FLOW CW AIRFOILS 


At sufficiently small airfoil Eeynolds numbers a region of 
separated flow is often found that has for its forward boundary the 
laminar separation point and for its roarrard boundary a tui-oiuent 
boundary layer (references 36 to 40 ) - The usual p.laoes- of -occur- 
ronce of the region of separated flow are near the leading edge of 
airfoils at angles of attack and behind the pressure 

point on airfoils wliich have extensive Icuninar boundary layers. 

Only q.ualit.ative information is avaijnble cQncerrd.iig the extent 
.of the region of separated florr. The flow can reattach Itself to 
the surface as a turbulent boundary layer at scaao smell dis’cance 
behind the point at wlilch the laminar botiadary layov leaves the 
surface. The extent of the region of separated frow lying between 
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the laMnar and turbulent bbxmdary layers has betn obwOVTed. to 
decrease -'With an increase iii Beynolds number (reference 40) . JSo ■ 
case is Joicwn for w*ich the separated turbulent' boundary layer iiao 
rejoined the surface ;ln free flow. 

The formation of turbulent Telocity f3.uctuations in the fluid 
layers that are moring do'wnstream and that are Just above the reg:’-'''n 
of separated flow is believed to be impor'tant in the foimation of a 
turbulent boundary :3eyer behind the region of separated flow. 


HJOFILE-ISBAG COMPOTATION 


The method of profile-drag coniputaticai is based on the m'.anentum 
theorem which may be. stated as follows: In a steady flow without 

body forces the net external force acting in a partlculaif direction 
on the surface bounding a fixed region of fluid is equal in magnitude 
to the difference between the time rate of outflow and time rate of 
inflow of momentum in 'the direction under consideration and has the 
sense in which the momentum decreases. 

The profile dz*ag can therefore be determi-ned if the difference 
between the time rate of outflow and time rate Inflow of momerit:mi 
in the direction of fllj^it can be determined for a region of fluid 
bounded by the bod^'' and a surface which has a shape and distance 
from the body so chosen that the ju*esBureB produce no resultant 
force on this surface 'In the line of flight direction . Squire and 
Young (reference Qk-) de'beimine "the profile drag by computing ■the 
momentum ■thiolaiess at the airfoil trailing edge from the von Karman 
momentum equationj then by the use of the von Karman momentum equa- 
tion together with an assumption concerrJ.ng the velocity profile 
across ■fche VQ’ke they compute -the momentum thickness very far behind 
the body. The profile-drag coefficient is kno^zn once the momentum 
thickness, and therefore tiie momentum defect, is known very far 
behind the body. 

The Squire and Yoimg method Is accurate to -V'Tithln a fe'\Y- percent 
if no regions of separated flow are present and if the tranoitlon 
point is known. It must bo eitrphasiged that the profile-di’feg coeffi- 
cient can be computed only when the transition point is kno-wii and 
that the profile-drag coefficient is sensitive to the position of 
the transition point. 
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PIEE FLOW 


The subject of flow tti pipes heo been studied ezhaustively and 
the knowledge of the subject may be found in many works ; for esnmple, 
references 1, 2, and 15 • The steady la m inar flow in a pipe can be 
computed directly from the eq,uations of motion. Included are the 
pressure-drop formulas and the velocity profile (reference 1, 
pp. 36-39). 

The knowledge of turbulent flow in pipes, like that of turbulent 
flow over bodies, is based on experiment. The pressure-drop formulas 
have been determined from flow experiments with smooth-wall and 
rou^-wall pipes (reference 1, pp. 135 “1^5; and references 2 and 15) * 
Universal velocity-distribution formulas have been determined for 
flow in smooth and rou^ pipes. 

The mixing-lengbh theories of Prandtl and von Earman make 
possible the con^jutation of the velocity dietributicax across the 
pipe ■vdien the surface friction is known. 


DIFBTOERS 


A diffuser is a duct having an Internal area that increases 
with distance downstream. Because of the increasing area, the flow 
velocity in the diffuser decreases with distance downstream and 
therefore the static presswe in the flow increases. 

A purely theoretical treatment of laminar flow in a two- 
dimensional diffuser is given in reference 2, which states that the 
results are of theoretical interest only. 

The literature for diffusers with turbulent flow is extensive; 
but because no theory of turbulent flow has been established for 
diffusers, the work of the various experimenters (for exampJ-e see 
reference 4l) has not resulted in the ability to predict the behavior 
of a given diffuser with a good degree of certainty. 


FREE -MIXING PROCESSES 


A free -mixing process is one in which solid boundaries play no 
part. Some examples of cases in which free-mixing processes occur 
are Jets, wakes, and regions in which parallel streams of different 
velocities meett The cases of turbulent-mixing processes, treated 
in references k2 to k5, are based upon an assumption relating the 
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local shearing etrass to the ch^acteristics of the local velocity 
profile. The eq,vetions and assmptione of these references are 

(1) The eq.uaticax of motion id.th zero static -pres sure gradient 

(far from, bodies) • • 

( 2 ) The equation of continuity 

f 

( 3 ) The Prandtl eq.uation for the shearing stress 

-,2 ou 
T » p2 ~ 
oy 

(k) The mixing length I proportional to the width of the 
mixing reglm 


5u 

cy 


(^0 


The cases treated in reference k6 use for the shearing stress the 
eq,uation 


T 



(5) 


instead of r 



The shaa:*pn-ese of the calculated 


velocity profiles at extreme values of u is eliminated by usljog 
equation (5) instead of equation (4) but no better understanding of 
the flow process results. 


Turbulent Jets and wahes and turbulent flows in the presemce of 
boundaries have also been analyzed by considering the velocity and 
pressure fluctuations in the flow. An introduction to this method 
of investigation is given in reference kj . 


EETDCT OF BOUWDAEY IAY1I?-0N P0TB3TOIAL-FL0W 
CHAEIACTERISTICS OF AIRFOILS 


Airfoil characteristics computed by potential-flow thooi’y are 
known to differ from the experimentally determined characteristics. 
The deviation from the potential-'flow c2aaract eristics appears not 
only in the existence of a profile drag but also in changes in the 
lift and pitching-moment characteristics. The changes are caused 
by the presence of the boundary layer and usually increase with 
Increasing boundary -layer thickness. 
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The hasis of the method of congjutation of refereme 48 is 
Taylor's theorem 'vjhich states that egual positive and negative 
amounts of vcsrticity are shed fi*om the airfoil trailing edge poi’ 
unit time when the lift is steady. The method of ccmputation of 
reference 48 in outline is approximately as follows; The potential- 
flow velocity distribution over the airfoil is found; a suitable 
fairing is made at the trailing edge if necessary to avoid a 
stagnation point. The boundary -layer thicknesses at; the trailing 
edge of the upper and lower surfaces are detejrmtned after choosing 
the transition point and the Eeynolds number. The velocities at 
the edge of the upper-surface and lower-sui*face boundary layers at 
the trailing edge are then cocrputed; vhen these, velocities are taioTO 
the pressure rise throu^ the boundary layer is estimated for both 
upper and lower surfaces. If the pressure at the trailing edge is 
not the same for both upper and lower sui’facoB, Taylor 'a tlaeorem 
is violated; the circulation is therefore adjusted and the procedure 
repeated until the pressures are e<jual. A.t present the oemputationB 
are too lengthy and, without empirical correction factors, are too 
inaccurate for routine iiae. 


Langley Ifemorial Aeronautical Laboratory 

National Advisory Conmlttee for Aeronautics 
Langley Field, Va., May 21, 1S>47 
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APP5MDIX 

MAIHIEHAriCAL DERIVATIOWS 
Prandtl Boundary-Layer Equatione 


The equations of motion in Cartesian coordinates for incom- 
pressitle flow are as folloT'®: 

The eqimtion of motion for ‘x-dii'ecti on, is 



The equation of motion for y~diroction is 



The equation of continuity is 


_ 

'dx by 


Eor the boundary layer, the equations become the Prandtl boundary 
layer equations; thus 


and 


P 



^ + i,2!s 
Sx 


(Al) 


0 = ^ 

du ^ 5v 

By 


(A2) 


0 
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Equations (Al) and (A2) are valid ^on the preasure is cosistanb 
across the thictaess of the houndary layer, when the ratio of tlie 
boundary-layer thiclaiess to the curvature of the surface is 
negligible, and -vdien all viscous terms iitvolving either v or 
derivatives with respect to x are negligible. Eor flows about 
airfoils at normal an^ea of attack and about plates at zero angle of 
attack these conditions are accurate over most of the length of 
the surface. Begions ■i&ere the ccaiditions may not be accxirate are 
in the vicinity of stagnation points and in the vicinity of the 
separation point. 


Derivation of the von Saroan Mcmientum Equation 



Continuity 

The mass leaving box abed equals the mass entering box abed. 
The mass entering per unit time throu^ ab is equal to 




pu dy 


The mass leaving per unit time through cd is equal to 


08 


08 


pu dy + As 


yo 


dx 


pu dy 
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Therefor©, .^e . Mass ©nterlfi^ wiit iime throuiflih be le ©q,iaal to 


H 5 




■ pu dy + ^ix — 
dx 

I/O ^ 00 


pu ay 



r »6 


pu dy) a Ax 


dx 


pu dy 


00 


Appllcsatlon of Mauentian Theorom 

The net change in moraanttttn. per ur^t time in the x-direction 
equals the net force acting to ri^t on box in the x-direction. 
The x-monentum entering per unit time throu^ ab is given by 

IS 

A 

pu*^ dy 

Jo 

The x-momentum entering per unit tine throu^i be is given by 


IT Ax — 
dx 



pu dy 


The x-momentum leaving through cd i'S given by 





pu^ dy + Ac ~ 
dx 


2 

pu dy 


Oo- 


Net Excess of x-Momentum Leaving to Ei^t 

over That Entering from Left 

The net excess of x-momentum leaving to ri^t over that 
entering from left is given by 


P 5 


pu^ dy +- Ax 

dx 


ns 


P 6 


PU^ dy 


Jo 


06 


OO 


pu^ dy - U Ax ~ 
dx 


pu dy 


UO 



KA.CA TN No, 1384 


27 


Net Forces Acting to Ei^t on Box 
The pressure acting to the ri^t on the box in given by 


pS + p — sin a, 
cos a 


The pressure acting to the left on the box is given by 


p8 + ^ /ix 
dx 


The sldn friction to the left on the box is given by 


Ax 


Therefore, the net force acting to the right on the box is given by 


dp6 

p8 + p Ax tan a " p8 - Ax - t Ax 

dx ° 


Then, since the not change in zaojasntuza in the X”diroction 
eq.uale net fox’ce acting to right on box, 


05 


Ax 


,dx 


ns 


tlo 


P d 

pir dy - U — 
dx 


pu dy 


Uo 


=j p Ax tan a-8-^Ax-p”Ax-T Ax 
dx dx ° 


dS 


There 


tan a = 


d& 

dx 
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Thon 


Ax 


Ldx 


f)S 


ns 


uo 


pu^ dy - U 

dx 


pu dy = 


ao 


\ dx o / 


or 


ns 


6^ 

dx 


Ps 


0 


pu^ dy - U ^ 
dx 


pu dy = 


do 


S - f 
dx 


(A3) 


After the definitions 


ns 


u^e 


u(U - u) dy 


do 


and 


ns 


US* = 


(U - u) dy 


0 


have teen used, eq.uation (A3) can he vritten as 


p 6 * ^ + T (Alt) 

dx dx ° 

If the eqiiation of motion for inviscid flow that is true outeido 

the boundary layer, - pU — , is .used and the term — is 

dx dx dx 

split into two terms, equation (A4) can he vcitten as 

dx \ U dx/ pxj2 
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The von Barman equation is often written in' this form. ■ Bj using 
the Bernoulli eq.'^tlon to relate the ntatic end d^nSamifC' pressure ■ 

outside the houndary layer, - ■§3'^ oq.ua tion (Ak) may also he 

dac djc 


written as 


^ + S + 2 £ ^ ^ 
dx 2 q. dx 2q. 


Derivation of the Esjiression fosp the Velocity Profile 
in the Bohlhaunen Method 

Let the houndary-layer velocity prorlle he given hy a fourth 
degree polynomial 


P ' *3 k 

u n a;y- + hy .+ cy-' + dy 

The values of the coefficients are determined from the farm 
taken hy the eq.uation of motion at the ■wall and from the assiuiption 
that all viscous effects are confined to a thin boundary layer. 

At the wall, u = v = 0, and so the eq.ua ti. on of motion 



Because all viscoxia effects are assumed to he absent outside 
the boundary layer, t = — = 0 outside the boundary layer. There - 

^ ' 2 
fore, at y = S with t = ^ it follows that ~ = 0. 

Sy dy 5 y 2 

Erom the definition of 5 it also follows that u = U. 
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By using thee© fovcc conditions, four simultaneous linear equa- 
tions ar© obtained fivan tdiich a, "b, c, and d are evaluated. The 
result is 


U(12 + X) 
65 



c 


U(4 - X) 



U(6 “ X) 
65^ 


■where the velocity profile shape parameter, equals 
expression for •tlx© velocity profile thoreforo beccmos ■ 


o dx 


The 


^ ^ ^ 4 - X /y\^ , 6 - X 

u ■ 6 W 5 \8/ 2 \sy 6 \B/ 


and the surface shearing stress 



is 


^o " 




U 12 + X 

a 6 
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I2 = i- + X 

2<1 Eg 6 


JrjL order to shov that X « - 


ixse tlxe definition of 


X o ^ S 
u &x 


It follows that 


_ i_ ^ :* ^ dp 5 dp 

■0 (SU dz pu2 dx ~ ® 2q. dx 


Then, from 


I2. - h, ^ ^ 

^ 9 . K5 6 


^ 2 q. 

E- « ~ 

O •n 

o 
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aM so 


X « - ^ ^ 

2<i dx 

or 


X « 


5 


dx 


T 


O 


To show that X is Independent of 


- (T) Jk M 

o dx \oJ dx/o u 



hut (see p. 13) 


therefore, 



X 



dx/c 


C onsequently \iben the pressure dlstrlhution does not change with 
X is independent of 



WACA TN No. 1384 


33 


EiS’SRENCES 


1. Prandtl, L.: The Mechanics of Viscous Fluids. Vol, III of 

Aerodynanic Theory, div. G, W.- F. Durand, ed. , Julius , ■ 

Springer (Berlin), 1935* 

2 . Fluid Motion Panel of the Aeronautical Research Ccianittee end 

Others; Modern Dovelqpiaents in Fluid DynaTnics . Vole. I 
and II, S. Goldstein, ed., Eie Clarendon Press (Oxford), 193S» 

3- von Eiarman, Th.; On Laminar and Turtulcait Friction., NACA TM 
No. 1092 , 1946. 

4. Fallcner, V. M.: A Further Investigation of Solutions of the 

Boundary Layer Equations. E. &, M. No. l884, British A. E.C. . 

1937 . 

5 . Falkner, V. M. : Siir^lified Calculation of the Laminar Boundary 

Layer. E. & M» No. I 895 , British A. E.C., 1941. 

6 . Ho'warth, L .5 On the Solution of the Laminar Boundary Layei- 

Eq.uations. Proc. Roy. Soo. (London), ser. A,> vol. l64, 
no. 919 , Feb. 18 , 193^^^ PP- ’547-579. . ■ 

7* von Karman , Th», and Millikan, C. B.: On the Theory of Laminar 

Boundary Layers Involving Separation. NACA. Rep. No. -504, 1934. 

- 8 . Ton Doenhoff, Albert E.: A Method, of Rapidly Estimating the 

Position of the. Laminar Separation Point. NACA TN No. 67 1, 

1938 . 

9 . Gortler, H.; Further Development of a Boundary Layer ■ Prof ile 
for a Given Pressure Disti’ibution. Jour. R.A.S., ’inl. XLV, 
no. 362 , Feb. 1941, pp. 35~50. 

10 . Wada, K. : Theory of Laminar Boundary Layer. Rep. No. I 96 

(vol. Sy, 16 ), Aero. Res. Inst., Tolyo Imperial Uni v.. Sept. 

1940 . • 

11. Loytzansliy, L. G.; Approximate Jfethod fcu’ Calculating the Laminar 

Bou n da r y Layer on the Airfqll. Ccmptes Eendus Acad. Scl. USSR, 
vol. 2SXV, no. 8 , 1942, pp, 227 .- 236 . . 

12. Eochin, N. E., and LqytzansJsy, L. G.! An Approximate Method of - 

Calcula-ting the Laminar Boundary Layer. CoraDtee Eendus Acad. 
Sci. US®, vol. S2ZV1, no. 9, 1942, pp. 262-266. 



NACA TN No. 1384 


34 


13* Falisner, V. M»J A New Law for Cfaiculating Drag. The Beeistance 
of a Smooth Flat Plate ^■dth Purbulent Boundary Layer. Air- 
' , craft Engineering, vol. XV, no. I 69 , March 1943^ pp« 65“69. 

14. Schultz “Grunow, F.: Now Fricticxnal Beeiatanco Law for Smooth 

Plates, NASA 2M No. 986, 1941. 

15. BaKhmeteff, Boris A.: The Mechanics of Turbulent Flow. 

Princeton tfaly. Press, 1936. 

16. Fage, A., and Falknar, V. M.: An Es^erimontal Determination of 

t^e Itttensity of Friction oii the Surface of an Aerofoil i 
B. & M. No. 1315, British A. E.C., I93I. 

17* Tcm Doenhoff, - Albert 1., and. Tetervin, Nealr Determination of 
General Eolations- for -the Behavior of Turbulent Boundary 
Layers. HACA ACE No. 3G13, 1943- 

18. Grtisch^.tz, E.: Die turbiileh-bb Eeibungsschicht in ebener 
Stromung bei Damckabfall und Druclcan.Btieg, Ihg. -Archlv, 

Bd, II, Heft 3, Sept. 1931, pp. 321-346. 

19* Fedioevsky, E.J Turbulent Boundary layer of an Airfoil. 

NACA TM No. 8B2, 1937. 

20. Kehl, A.J Investigations- on Convergent and Divergent Turbulent 

Boundary Layers. B.T.P. Translation No. 2035.» British’ 
Ministry of Aircraft Production. (From Ing. -Archiv, Bd. 13^ 
Heft 5,. 1943, pp. 293-329.) 

21. Kalikham, L. E.; A New Ifethod for Calculating the Turbulent 

Boundary Layer and Determining the Separation Point. Comptes 
Eendus Aoad. Sci. USSR, vol. XXXVIH, no. 5-6,' 1943, ' ' 
pp. 165-169. 

22. Schmidbauer, Hans; Behavior of Turbulent Boundoi-*y Layers on 

Curved Convex Walls. NACA TM No. 791, I936, • 

23" Peters, H.; On the Separation of Turbulent Boundary Layers. 
Jour. Aero. Sci., vol. 3, no. 1, Sept. 1935, pp* 7“12. 

24. Sq,iUre, H. B., and Young, A* D.; The Calculation of -fche Profile 

Drag of Aerofoils, E. & M. No. 1838, British A.R.C,'1938. 

25. Tetervin, Neal; A i4e-f^gd for the Bapid Estimation of Turbulent 

Boundary -Layer Thicknesses far Caloulatiiag Profile Drag. 

NACA ACE No. L4G14, 1944. 



NACA TN No. 1384 


35 


26 . JacolDe, E. N., and von Doenhoff, A. E.;- Formulas fear Use in 

Boundary-Layer Calculations cn Lov->3)rag Wings ► NACA ACR, 

Aug. 1941 . 

27 . Schlichting, H.: Esperimental Investigation of the Prohlem of 

Surface Eoughness. NACA TM Ko» 623 , 1937* 

28 . Colahrook, C. F-, and White, C. M.: Exjieriments -with Fluid 

Friction in Eou^enad Pipes . Proc. Eoy. Soc. (London), 
ser. A, vol. 16 I, Aug. 1937> PP* 387 - 3 S 1 . 

29 * Iheodoi'sen, Theodore, and Eogier, Arthur: Experiments on Drag 

of Eevolving Dlslca, Cylinders, and Streamline Eods at Hi^ 
Speeds. NACA ACE No. L4F16, 1944. 

30 . von Karmln, Th.s Ttarh?J.ence and Skin Friction. Jour. Aero. 

Sci, vol. 1, no. 1, Jan. 1934, pp. 1-20. 

31 . Jones, B. Melvilli Fli^t Experiments on the Boundary Layer. 

Jour. Aero. Sci, vol. 5^ no. 3, Jan. 1938^ PP* 81 - 94 . 

32 . Becker, John V.: Boundary-Ia.yer Transition on tl:e W.A.O.A. 0012 

and 23012 Airfoils in the 8-Foot Hi^-Speed 'Wind Tunnel. 

NACA ACE, Jan. 1940. 

33 » Liepmann, H. W.; Investigation of Boundary Layer Trar.sitlon on 
Concave Walls. NACA ACE No. 4J28, 1945 . . . ' . 

34 . Lin, C. C.: On the Stability of T-sTO-Dimensional Parallel Flows.’ 

Paxi; I. Qxiarterly Appl. Math., vol. Ill, no. 2, July 1945, 
pp. 117 - 142 ; Part II, vol. Ill, no. 3, Oct. 1945, pp. 218-234; 
and Part III, vol. Ill, no. 4, Jan. 1946, pp. 277"301. 

35* Schu’bauer, G. B., and Skramstad, H. K.; Laminar -Boundary -Layer 
Oscillations and Trai^ition on a Flat Plate. NACA ACE, 

April 1943 . 

36 . Schlichting, H., and Ulrich, A.; Zur Berechnung dos Umschlages 

laminar turbulent. Jahrbf 1942 der deutschen Luftfahrtforschung, 
E. Oldenbourg (Munich), pp. I 8 - I 35 . 

37 • Pretsch, J . : Die Stabllitat der Laminarstrmung bei Druckgefalle 

und Druckanstleg. Fersohungsbericht Nr. 1343, Deutsche 
Luftfahrtforschung (Gottingen), 1941. 

38 . Jones, B. Melvlll: Stalling. Jour. E.A.S., vol. jcmrill, 

no. 285 , Sept. 1934 , pp. 753 - 770 . 



36 


WACA TN No.. 1384 


* 


v/ 39* Schubauer, G. -B.: Air Flow in the • Boundary Laj-er of an Elliptic 

Cylinder. NACA Eep. No. 652/ 1939* 

40. von Doenhoff, Albert E.^ and Tetervin, Heal: Invoetigation of 

the Variation of Lift Coefficient AvLth Rejhiolda Hiaiber at a 
Moderate Angle of Attack on a Low-Lrag AirfeiJ.. NACA CB, 

Nov. 1942. 

41. Patterson, G. H.; Modern Diffuser Design. Aircraft Engineering, 

vol. X, no. 115, Sept. I938, pp. -267-273. 

« « 

42. Schlichting, H.I Uber das ebone Wind-schatten-problcan. Ing. "Ar'-’V; j v, 

Bd. I, Heft 5, Dec. 1930, pp.- 533“571. 

43. Tollmien, Walter: Calculation of Turbulent Espansion Processes. 

WACA TM Ho., 1085, 1945, ■' 

44. Forthmann, E.: Turbui.ent Jet Ispansion. WACA TM Ho. 1936. 

45. Kuethe, Arnold M.; InvestigatlcaoB of the Turbu-leht Mixing 

Regi-ons Formed by Jets. Jour. Appl. Mech., vol. 2, no. 3, 

Sept, 1935, pp. A-&T A-95. 

46. Gortler, H.: A Hew Approximation Method for the Nuxaerical 

Evaluation of Free Turbulence. Problems . R.T.P. Translation 
Wo . 2234, British 14Lnistry of Aircraft Production. ■ (From 
Z.f.a.M.M., Bd. 22, Heft 5, Oct. 1Q42, ppV 2it0-254 . ) ' - 

47. Dryden, Hu^^ L.: A Review of the S-fcatistical Theory of Turbulence. 

Quarterly Appl. Math., vol. I, no. 1, April 1943> pp.-7-42. 

48. Preston, J. H.: The Approximate Calculation of the Lift of 

Symmetrical Aerofoils taldiig Account of ’the Boundary Layer, 
with Application to Control Probleias. R. &"M. Wo. I996, 

British A.R.C., May 20, 1943. . ‘ ’ 








NACA TN No. 1384 


Figs. 5,6 



NATIONAL ADVISORY 
COMMITTEE FOfi AERONAUTICS 


Figure 5«- Boundary-layer thickness on upper surface of airfoil. 
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Figxire 6.- Flat-plate turbulent and laminar velocity profiles 
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Figure ?•- Effect of adverse pressure gradient on velocity profile 
of turbulent boundary layer* 
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